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ABSTRACT 

Analytic expressions are given for the r decays into three charged leptons, r — > ItlvrV^ 
where the I are combinations of electrons and muons, and for the radiative decays r — > ^Iv^n. 
The branching ratios are sensitive functions of whether the final state I are muons or electrons. 



The branching ratio for the process fi + — > e + e + e~u e V^ was first measured more than 
thirty years ago 1-4 . The measurements are in good agreement with the standard model 
theoretical predictions 5 . More recently, the primary interest in this process is as a background 
for processes that would signal new physics beyond the standard model such as /i + — > e + e + e~. 

The study of r decays is rapidly reaching the point where similar studies could be done. At 
CLEO, for example, enough rs' may have already been produced to allow measurement of the 
branching ratios for the four leptonic final states, e + e + e~u e V r , n + e + e~u^U T , e + fj, + \i~vJj t 
and fi + fi + \i~v^ T \ or to place limits on the branching ratios for the corresponding lepton 
number violating processes which do not involve neutrinos in the final states 6 . The purpose 
of this note is to give the standard model branching ratios for £ + £ + £~z/^I7 T , where I = e or 
/i; and to present an analytic expression for the square of the matrix element which can be 
easily incorporated into Monte Carlo simulations. 

The processes proceed through the Feynman diagrams given in Figure 1 where the last 
two diagrams are needed only for the cases of identical particles in the final state. Each term 
in the square of the matrix element is of the form 

Tr /(l_ 75 ) 7 .p 57 «(l_ 75 )Ax Tr la (1 - 75 ) 7 -p ilp (1-75)5 

where A and B are strings of gamma matrices and p 4 and p 5 are the momenta of the neutrinos. 
If these traces are done directly the result involves many terms. However, the result can be 
made shorter by first introducing an identity matrix on either side of A and of B, and Fierz 
rearranging both. The expression then becomes 

Tr A la (1 + 75 ) Tr 7/3 (1 - 75) 7 • p 5 7a (1 - 75) 1° 
xTr £7,(1 + 75) Tr 7 a (1 - 7s ) 7 ■ p 4 ^ (1-75)7" 
~ p p 5 pi Tr A la (1 + 75) Tr B lp (1 + 7s ) 

and the squared matrix element is greatly simplified. Further, since this is the only dependence 
on P4 and p$, the phase space integrals of the neutrinos can be easily done. 

The total rate for the process t + (P) — > ^ + (pi) + e + (p 2 ) + e~(p 3 ) + V T + is given by 



r = 



G 
72 



where = P M — pi — p% — p%, P 2 = M 2 , p\ = m 2 p\ = m 2 e = p§, and s = (p 2 + P3) 2 . 
The functions are what remain after analytically computing the integrals over the neutrino 
momenta. The rest of the squared matrix element, T, is given in Appendix I. Expressions 
for the rate to the other final states can be obtained by interchanging m M and m e or by 
setting them equal. In addition, for the identical lepton cases, the expression for \T has 
to be generalized as explained in the appendix. The rate for /i + — > e + e + e~ Vep ^ can also be 
obtained from (1) by setting to m e and then M equal to m^. 
The branching ratios are given by the following table 





-> e + e + e 




3.60 ±0.02 x 10- 5 


r+ - 




~e + e 


Z/ e 17 r 


4.15 ±0.06 x 10~ 5 


r+ - 


-> fi 






1.97 ±0.02 x 10~ 5 


r+ - 


■* 


>v 


~u e V T 


1.257 ±0.003 x 10 


r+ - 


-> i£ 






1.190 ±0.002 x 10 



where the errors are due to inaccuracies in the numerical integration. These numbers assume 
the branching ratio for r — > eu r u e is 0.177. 

For the tau cases where the pair produced are electrons, the branching ratio is of the same 
order as the muon decay. If the pair produced are muons, however, the branching ratio is 
two orders of magnitude smaller. This seems to be just the effect of the 1/s 2 in Eq. (1); for 
electrons s can be a lot smaller than for muons. (s > 4m 2 where m is the electron or muon 
mass.) This means that any experimental cut on the minimum value of invariant mass of the 
produced pair will drastically reduce the number of events. This also means that the cross 
terms between the first and second diagrams of Fig. 1 and the third and fourth diagrams 
make only a small contribution to the rate because these cross terms go as 1/ss', where 
s ' — (pi ±P3) 2 ) an d s an d s' cannot be small simultaneously. These observations were verified 
by our numerical results. 

An analysis of the angular distributions of the final state leptons reveals strong correlations. 
For the case with three electrons in the final state we have plotted, in Fig. 2a, the distributions 
in the angles between the two like-sign and the opposite-sign electrons. In both cases the 
distributions peak when the two momenta are parallel. Similar results are obtained for the 



2 



final state n + e + e~u f J7 T , but the correlations disappear for the final states e + u + [i~ v e V T and 
/i + /i + yTv v p T . Again these results are simply a consequence of the 1/s 2 term dominating the 
dynamics. In Fig. 2b we have plotted the distributions in the total energy E c of the final 
state charged leptons. For comparison we also display the same distribution for just the phase 
space without any dynamics. The dynamics shift the E c distribution to lower values. The 
shift to lower values increases with decreasing mass of the final state leptons. 

For completeness we also include an expression for the radiative decay r(P) — > t{$>-i) + 
j(k) + v T + vg for i = e or /i. With the neutrinos again integrated out the rate is given by 

F = ^^/^^0(^)e(^)T 7 (2) 
(2vr) 2 3Mi 2w 2£ 2 1 M M K J 

where = P M — k^ — -p^ and T 7 is given in Appendix II. The photon must be visible and 

thus the photon energy has some nonzero lower limit uj§. This form would seem to be useful 

if further experimental cuts need to be applied to the photon or charged lepton. If, however, 

no further cuts are required, then the integrals can be done analytically and the branching 

ratio to r — > i + v T + vg is given by 7 
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3tt 







{(<■«-- 
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12 J 



6 in 



Do 



+ 3[L(l)-%)]--[6 + (l-|/ ) 3 



6(l-y )-(l-yo) 4 

1 - y )£n(l - y ) 



+ 1 (1 - 2/ )(125 + 45y - 33y 2 + 7y 3 )} , 



(3) 



where y = w /wMax and L is a Spence function defined with an overall sign such that 
L(l) = — 7r 2 /6. A graph of the branching ratio given by this equation is shown in Fig. 3 for 
the three processes /i^e + 7 + z/ + 77, r^e + 7 + z/ + 77, and r — > /x + 7 + v + V. 

In conclusion we have found that with about one million rs' the standard model branching 
ratio for r decays to the final states e + e + e~u e U T and n + e + e~u^U T may be measurable. The 
distributions in the angles between the final state leptons are strongly peaked at zero and 
stringent angular separation cuts can drastically reduce the rates. 
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APPENDIX I 

For the decay t + (P) 



/i + G»i) + e + {p 2 ) + e (p 3 ) + i/ T + ^ define 





— r Pi P2 


e 


— \P2 P3) 




— {Pi T P3j 




= P1-P2 


■u 


= Q-P 


Sl 


= P-Pl 


«2 


= P-P2 


S3 


= P-P3 


«1 


= Q-pi 


U 2 


= Q-P2 


u 3 


= Q-Ps 



Obviously there are relations among these invariants but the following analytic expression is 
simpler if we use an overcomplete set. In terms of these the expression for T in Eq. (1) is 
given by 

T^j— : r [ - {Q 2 si + 2uui)(m 2 e + s 2 + s 3 ) 

s — 2{s 2 + s 3 J 

+ - Q 2 (2s 2 t + s 3 (s' -m 2 e - m 2 )) + 2wi(w 2 S2 + u 3 s 3 ) 



T = 



[s-2(s 2 + s 3 )]< 



-Q 2 (2 Sl - 2t - s' + m\ + mj) + 2u x {u - u 2 - u 3 ) 



x 



4 + 4 + -sM 2 + m 2 e (s 2 + s 3 ) 



+ 



1 



s + s' + 2t — ml — 



(S Tfl \ 



+ (Q 2 s 3 + 2mm 3 ) - m 2 . - mj) 



+ 



s + s' + 2t — m\ — m^ 



Q 2 (si + s 2 + s 3 ) + 2-u(«i + w 2 + M3) 



X 



-t 2 - i(s' - m\ - ml) 2 - i sml + \m 2 e (2t + s' - m\ - m 2 ) 
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+ s - 2(s 2 + S3) s + s ' + 2t - ml - 777,2 
x - s(s 2 + s 3 )(Q 2 ml + 2u 2 2 ) 

+ (Q 2 t + 2m!M 2 ) (ssi + s 3 (s' - 2t - s + m 2 - mj) + 2m 2 e s 2 ) 
+ X - (Q 2 (V -ml- ml) + 4uiu 3 ) 

x (ssi + s 2 (2t - s' - s + ml + 3t77 2 ) + 2m 2 e s 3 s ) 

+ 2 (Q 2 Si + 2u«i) (2tS 3 + S 2 (V - 777 2 - 7T7 2 ) - S (s X + 777 2 )) 

+ 2 (Q 2 777 2 + 2M 2 ) (s 3 (V - 777 2 - 777 2 ) - m^) 

+ 2 (g 2 777 2 + 2 U fj (2tS 2 - 777 2 Si) 

+ (q 2 ^ - Q 2 m\ + 2t7 2 77 3 ^ (2ssi - 4tt7 2 si - 4ts 3 - 2s 2 (s 1 - m 2 - mj)) 

+ (<5 2 S 2 + 2UU 2 ) ^(S 2 - S3) (s' - 777 2 - 777 2 ) + 2777 2 t - SSi + ^ (2777 
+ (g 2 S 3 + 2t777 3 ) (2t (s 3 - S 2 + 27T7 2 ) - S («i + t) + 777 2 (s' - 777 2 - 777 2 )) 
+ X - S (Q 2 M 2 + 277 2 ) (2t + S / - 777 2 - 777j) 

+ 6Q 2 [-ml (s 2 + s 3 ) (2t + s' - m 2 e - m 2 ) + ss 3 t + 2m 2 e s Sl + ^ss 2 (V - m 2 - m 2 ) 
If all final leptons are the same flavor set m 2 = m 2 and replace -^T by 

1 1 ™ 1 1 ^ / A 1 1 1 rr,, 

2^ + 2^ (Sl ^ s 2,«i^ 2 ,s^s) + - --T 

where T' is given by 

I I 



l-s) s'-ml 



s-2(s 2 + s 3 ) s'-2(s! + s 3 ) 

X 



(Q 2 t + 2t7i77 2 ) (4s 2 - 2m 2 M 2 ) 
+ [Q 2 - m 2 )j + 2 Wl u 3 ^ (-2s 2 s 3 + s (s 3 + M 2 ) - m 2 [2M 2 + s 2 + s 3 )) 
+ (Q 2 si + 2t777i) ^2s 2 s 3 - ss 3 - \-M 2 s + 27T7 2 (s 2 + S3) 
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+ (q 2 Q - m£) + 2m 2 m 3 ) (-2sis 3 + s's 3 + s'M 2 - m 2 (2M 2 + s x + s 3 )) 

+ (g 2 s 2 + 2ww 2 ) (2S1S3 - s's 3 - i AfV + 2m 2 ( Sl + s 3 ) 

+ (g 2 m 2 + 2m 2 ) (-2t (s 3 + M 2 ) - m 2 (si + s 2 )) 

+ (g 2 s 3 + 2ww 3 ) (M 2 + 4s 3 ) + m 2 ^2s 2 + 2s x + t + -s + ^s' - M 

+ (<Q 2 M 2 + 2m 2 ) (-2ts 3 + m 2 (2s 3 - 2* - s - s' + 2m 2 ) ) 

+ 6Q 2 (-2ts 2 + m 2 (-2s 2 s 3 - 2 Sl s 3 - 2 Sl s 2 + ^s' (s 2 + s 3 ) + ( Sl + s 3 ) 
+ M 2 (2t + s' + s-2m 2 ))) 



5 m 2 



[s + s' + 2i - 2m2] : 



x 



(<5 2 (si + s 2 ) + 2u(ui + m 2 )) (-2t(s + s') + 2m 2 (s + s' - 2m 2 )) 
+4 (g 2 s 3 + 2mm 3 ) i (f - 2m 2 ) 



s + s' + 2t - 2m 2 e s - 2(s 2 + s 3 ) 

x (Q 2 m 2 + 2m 2 ) ss 3 + (Q 2 m 2 + 2m 2 ) (-s's 3 + m 2 ( Si + s 3 )) 
+ (g 2 m 2 + 2w 2 ) (-2ts 2 + m 2 ( Sl + s 2 )) 
- (g 2 M 2 + 2M 2 )st 

+ (QH + 2 Mi m 2 ) (s 3 (s - s') - m 2 (s 2 + 3s 3 )) 

+ ^g 2 ^ - mlj + 2miM 3 ^| (s's 2 - ssi + m 2 (s 2 - s 3 )) 

+ (g 2 si + 2mmi) (-2ts 3 - s's 2 + ssi + m 2 (s + 4s 3 - 2s 2 )) 
+ (q 2 Q - m 2 ^ + 2w 2 w 3 ^ (2ts 3 + s's 2 - ssi + m 2 (2si - s 2 - s 3 

+ (g 2 s 2 + 2mm 2 ) ^ssi - s's 2 + m 2 ^2s 3 - - ^s' - t + m l^j 

+ (g 2 s 3 + 2ww 3 ) (t(s + 2s 2 ) + m 2 (2s 3 - 2s 2 + ^s-^s' -t + m 2 )) 
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+ 6 Q 2 (-sts 3 + m 2 e (t (s 2 + s 3 ) + ^s' (s 2 + s 3 ) + (s 2 - s 3 - 2si)) 



- m e (s 2 + s 3 
1 1 



s + s> + 2t - 2m 2 e s' - 2(si + s 3 ) 
(Q 2 m 2 e + 2ufj (-ss 3 + m 2 (s 2 + s 3 )) + (Q 2 m\ + 2uf) s's 3 

+ (Q 2 m 2 + 2uf) (-2tsi + m 2 (si + s 2 )) - (Q 2 M 2 + 2m 2 ) s'f 
(g 2 t + 2 Mi m 2 ) (s 3 (s' - s) - m 2 ( Sl + 3s 3 )) 



+ 



+ Q 2 



— - m 2 J + 2-uim 3 J (2ts 3 - s's 2 + ssi + m 2 (2s 2 - si - s 3 

+ (ffsi + 2uu 1 s j (^s's 2 - ssi + m 2 e (2s 3 - \s - ^s' -t + m 2 )^ 

+ (Q 2 Q - m^j + 2u 2 u 3 ^j (ssi - s's 2 + m\ {s x - s 3 )) 

+ (Q 2 s 2 + 2ww 2 ) (-2ts 3 + s's 2 - ssi + m\ (s' + 4s 3 - 2si)) 

* (V + 2si) + m 2 ( - t + 2s 3 - 2si + -s' - -s + m 2 ) 

+ 6Q 2 ( - s'ts 3 + m 2 (t( Sl + s 3 ) + ^s'(si - 2s 2 - s 3 ) + ^s(si + s 3 )) 

+ (-«! - s 3 ) 



APPENDIX II 

The decay t(P) — > £(p 2 ) + 7(fc) + ^ + 77 is given in the text in terms of T 7 . If the energy 
of the charged lepton, the energy of the photon, the mass of r, and the mass of t are denoted 
by E 2 , cu, M, and m, and k ■ p 2 is named u, then T 7 is given by 



T 7 = 4m 2 



1 1 

+ 



+ 



+u 

1 

u 



- ( -5M + 3— - 12P 2 ] + At f 3M 2 - 8MP 2 + 3m 2 ) - 4 

- 4w 2 M 2 + w (3M 3 - 12M 2 P 2 - 5Mm 2 ) 
+ - (6M 3 E 2 - 4M 2 m 2 P 2 - 8M 2 E% + 6Mm 2 P 2 ) 
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+6M 3 E 2 + 3M 2 m 2 - 16M 2 E 2 - 2Mm 2 E 2 + 3m 4 



+ 



AM 2 m 2 uj 2 + uj (-3M 3 m 2 + 8M 2 m 2 E 2 - 3Mm 4 ) 



-3M 3 m 2 E 2 + 2M 2 m 4 + AM 2 m 2 E 2 - 3Mm 4 E 2 

+ 4Mu + 2M 2 + 16ME 2 + 2m 2 
1 

+ - 



3M 3 + 2M 2 E 2 - 3Mm 2 + 16ME 2 - 6m 2 E 2 



+ 



0J' 



-3M A E 2 + 2M 2 m 2 + AM' Z E Z 2 - 3Mm'E 2 



2 77>2 
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FIGURE CAPTIONS 



1. The Feynman diagrams for the decays t + — > £ + £ + £ vtv T where the t are electrons or 
muons. 

2. (a) The distribution in z 12 — P\ • Pi for the two-like sign electrons (dashed line) and 
in Z13 = pi ■ p3 for the two-opposite sign leptons (solid line). For comparison the 
corresponding distribution for the phase space is also plotted (dash-dot line). 

(b) The distributions in the total energy E c = (Ei + E 2 + E 3 )/m T (solid line), and the 
corresponding distribution for the phase space (dashed line). 

3. The branching ratio, R, for /i — > e + 7 + u + v, r — ■> e + 7 + 1/ + u and r — > fj, + 7 + z/ + v as 
a function of the minimum detected photon energy, u . y = uj /u Msk . Note that R is 
defined relative to the three body process fi — ■> e + z/ + I7, r — > e + z/ + 77or r — > yU + z/ + I7 
and thus, for the r decays, is not the absolute branching ratio. 
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This figure "figl-l.png" is available in "png" format from: 



http://arXiv.org/ps/hep-ph/9402262vl 



This figure "figl-2.png" is available in "png" format from: 



http://arXiv.org/ps/hep-ph/9402262vl 



This figure "figl-3.png" is available in "png" format from: 



http://arXiv.org/ps/hep-ph/9402262vl 
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